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Domain walls in the spin S quantum Ising chain at zero temperature are consid- 
ered. Quantum effects prevent the walls from being sharp and delocalize the exact 
eigenstates. The absence of a characteristic length of the domain walls is reflected 
in the finite-size scaling of the spectrum of the Hamiltonian. 



1 Introduction 

Recently, there has been a growing interest in the study of interfaces with 
nontrivial geometry. Such interfaces may arise for example in domain walls 
in random magnets, fluid invasion in porous media, spreading on heteroge- 
neous surfaces, biological membranes and vesicles or epitaxial growthEJ When 
do quantum effects play a significant role in such problems ? For static prop- 
erties in non-random systems such as an spin S antiferromagnet with nearest 
neighbor exchange interaction J, the following regimes exist. Near the critical 
temperature Tjv ~ JS 2 , thermal fluctuations dominate. In the ordered phase, 
for temperatures T c /S <C T -C T c , quantum effects due to the finiteness of S 
are unimportant. For T < T c /S, however, one is in the quantum regime, where 
the quantum statistics of spin waves and their interactions lead to dependences 
on S and T not present in the classical (S — > 00) limit. 

Useful insight to the problem may be gained by considering a simple ID 
quantum system. We shall consider the domain walls at zero temperature of 
the following extension of the spin S Ising quantum chain 

H = "I? E 5 » - 2^ E S z {n)S z {n + 1) - ±L S y (n)S y (n + 1) (1) 

n— 1 71— n— 1 

where S x , y ,z(n) is a quantum spin S operator at site n, H is the transverse 
field, e describes the anisotropy in the spin space, J is the exchange coupling 
and L is the number of sites. We shall work with the boundary conditions 

S Z (0) = -S Z (L + 1)=S (2) 
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Table 1: Ground state phase of the quantum Hamiltonian W. The phase boundaries are for 

spin S = s. 



phase 


condition 


m 


G(i?) 


D 


e < ft- 1 





fl-V a exp(-i?/e) 


F 


ft - 1 < e < ft 2 / 4 




i?- 2 exp(-i?/C) 





ft 2 /4 < e 




iT 2 exp(-2i?/05R(Se aKi? ) 



The study of finite values of 5 rather than considering the classical limit S — > oo 
may be oLinterest, since several experimentally studied systems have a rather 
low valucBa of S: S = \ in CsCoCl 3 and CHAB, S = 1 in CsNiF 3 and S = § in 
tmmc. Emfiles of various physical quantities have been investigated by several 
authors H'BoB'H, Here we shall concentrate on the local magnetization profile 

m(r) = S- 1 (0\S z (n)\0) , r = n/(L + l) , n = 0, 1, . . . , L + 1 (3) 

where |0) is the ground state of TL. In the sequel, we shall often use the notation 
h = H/J. 

Depending on the parameters h and e, there are three distinct,ground state 
phases as indicated in table [l]. Their properties are distinguished!] through the 
magnetization m = ^ r m(r) and the connected two-point function G(R) = 
(S z (R)S z (0)) - (S z (R))(B z (0)). The correlation length £ and the wave vector 
K are for S = | knowifl functions of h and e. For higher spin S > ^, the 
F/O boundary is exactly known|3 while the D/F boundary (which is in the 
Ising universality class) has to be found numericallyEl Besides studying the 
spin dependence of the profile eq. (||) , we shall also ask how the profile changes 
between the two ordered phases F and O. 

Finally, the results presented here may have a bearing on certain non- 
equilibrium reaction-diffusinp processes, when the master equation is rewrit- 
ten in a Hamiltonian formE3 The resulting quantum Hamiltonians are almost 
identical to eq. (jl|) (at the O/F boundary) and surface fields model particle 
injection and extractionEJ Reactions of a single particle species are described 
by spin |, where the free fermion condition means that two particles react with 
infinite rate on encounter. That condition is apparently realised for the exciton 
kinetics in tmmcEj Ground state order parameter profiles as considered here 
should be analogous to particle density profiles in the steady state. 

In section 2, we review the results for the order parameter profile. In 
distinction to the classical case, quantum effects lead to the derealization of 
the energy eigenstates and thus to wide profiles. In section 3, related spectral 
properties of Ji will be described. 
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Figure 1: Magnetization profiles m(r) for spin S = h (boxes), 5 = 1 (full circles) and S = S 
(open circles). The curve gives the profile of eq. (^) which is correct to leading order in h. 

2 Order parameter profile for S = 1/2, 1 and 3/2 

We now display the results for the profile eq. (||). On finite lattices, the cal- 
culation of the profile from the lowest eigenstate of TL through the Lanczos 
algorithm is completely standardta We begin with the profiles for the phase F, 
which is shown in Fig. |l| for S = |, 1 and |. We observe the followingQ First, 
data for different lattice sizes (we used lattices up to L — 16) nicely collapse 
onto each other which implies that for the parameters considered the lattices 
are large enough to faithfully represent the L — > oo limit. Second, the profiles 
are quite wide and vary continuously with r. Third, the form of the profile 
is, at least for S not too large, quite independent of the value of S. These 
differences -grow larger when h approaches the critical value (if e = 0, h c = 1 
for S = ±P h c ~ 1.32587(1) for S = fl). Finally, the profile only depends 
weakly on h for h small enough. 

Finding a quantum profile with a width w ~ L of the order of the system 
size L is in contrast with the sharp profiles of width w ~ a, the lattice constant, 
in the classical S — > oo limit. This result can be explained as followscrQ The 
boundary conditions eq. (^) require that at least one domain wall is present in 
the system. For h = (and e = 0) , the states of lowest energy of Tt are 

ITT •■■!!) IT ---n •■■!.) IU---U) (4) 
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for spin S = | and for L sites, there are L + 1 of these degenerate states. For 
5 = 1 at h = the (L + 1) + L = 2L + 1 states 

ITT •••TO IT ••• Ti ••• i) ITI •••II) 
ITT ••■TO |> IT - - - T o i - - - i) |T o | •••ll) (5) 

are degenerate. These states have a single sharp domain wall. Turning on h, 
the transverse field term acts like a hopping matrix between these states and 
the domain wall profile becomes broad as a result of the superposition of many 
sharp domain walls. Working in the subspace M. of the states eqs.. (f|,|[) with 
a single domain wall, the order parameter profile is found to be&til 

m(r) = 1 — 2r H sin(27rr) (6) 

for both S = | and S = 1. As can be seen from Fig. [IJ for h small enough 
eq. (^|) is in excellent agreement with the data. 

Indeed, for spin S = \, the correspondence between the domain walls and 
the lowest excitations can be made precise&u. Namely, for the quantum chain 
eq. (|l|) with open boundary conditions one can define the domain wall creation 
operators, with a n = h (S x (n) + iS y (n)) 



D n = \ ( c «+i +c n +i) + ^ (ct -c„) , c„ = exp 



in a\att 



(T) 



which satisfy the anticommutation relations {D+ , D m } = <5„ !m , {D^,-D+} = 
0. and D„, respectively, create and annihilate a domain wall between the 
sites n and n + 1. Then, for open boundary conditions rather than the ones 
specified in eq. (|^) and for S 1 = ~, the quantum Hamiltonian becomes □ 

L-l . L-l 

W = J E D » D « + 2 eJ E (At - D+_ ! £>+ +1 + h.c.) 

n— 1 n— 1 



l L i 

- 2^ E ( D n-i Dn - D t-i D t + h.c.) - ^ J(£ - 1) (8) 



n=l 



where _Di = Z?o = ^(c^+ci+cj — c^). The Hamiltonian thus describes domain 
wall transfer between nearest neighbor and next-nearest neighbor sites as wall 
as domain wall pair creation and annihilation. This gives a nice interpretation 
to the well-knowro free fermion formulation. While in the F phase the next- 
to-nearest neighbor terms are unimportant, this is not so in the O phase. 
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Figure 2: Magnetization profile in the O phase for spin S = ^ and h = 0.01, e = 0.5 for 
several values of the system size L. 



Turning the profile in the 0_phase, different results are obtained from finite 
lattices for L even and L oddQ For L even, the lowest ground state remains 
non-degenerate and eq. (||) can be used. That is not so for L odd, where there 
are for h — and e ^ two degenerate ground states |0) , |0') and these states 
remain close together even for h ^= 0. Then the order parameter profile is 
found from the eigenvalues m± (r) of the matrix 



finite-size effects are much larger than in the F phase. Second, m(r) displays 
for L finite step- like behavior and it looks as if the system were built from hard 
objects each occupying two lattice sites. When h is increased, these composites 
soften until they melt at the O/F transition. Third, for L odd, the reflection 
symmetry is lost. These results can be explained □ in the same way as before 
through the hopping of domain walls, as described by the Hamiltonian in the 
form eq. (||). Calculations become particularly simple for h — and in the 
context of the subspace M. of the states eq. (0). We can then decompose 




(9) 



and m+(r) and m_ (r) = — m + (i — r) are related to each other through-spatial 
reflection. One of these is shown in Fig. 0. We observe the following^ First, 



M. = A4 e + M. into the subspaces with an even and odd number of inverted 
spins, respectively. Then, the effective Hamiltonian decomposes into two which 
are identical to the one found in the treatment of the F phase, but for a 
redefined effective system size. We have L e ^ — + for L odd (twice) and 
-^eff = L/2 + 1 and L/2 for L even, ft follows that in the L — > oo limit the 
width of the terraces decreases and the order parameter profile goes toward a 
smooth and reflection symmetric limit. The phenomenon seen in Fig. || is thus 
a finite-size effect. 

Finally, in the L — > oo limit, we observe (for spin S = =) numerically 
a simple relation of the profiles m(r; h, e) between the F and the O phases, 
suggesting thatQ 

m(r;/i eff (e),0)~m(r;0,e) (10) 

In the limit e -C I, the restriction to the subspace M. is accurate and yields 
/i e g-(e) = e. For e finite, we find that phenomenologically h e g(e) ~ sfe. 

With increasing values of S the mechanism of quantum derealization 
through hopping of the domain-walls becomes less and less important as the 
classical limit is approached W&B Again, this is best understood within the 
subspace M of states with a single domain wall. The quantum hopping of a 
domain wall becomes analogous to a tunneling process, since there is at least 
in the semiclassical limit S — > oo a finite energy barrier which prevents the 
free motion of the domain walls. This energy barrier vanishes for finite values 
of S in the thermodynamic limit. Starting from S = |, the states can be 
such arranged that the hopping matrix elements betwee±L.them become peri- 
odic functions, thus leading to a banded energy spectrumE_a The width of these 
bands decreases rapidly.™ with S, W ~ H2~ s . Finally, for £ = 00, quasiclas- 
sical approximations BB become accurate and yield a sharp domain wall (for 
h <C I) of the order of the lattice spacing. 

The effects of temperature and of Lohomogeneous transverse fields on the 
quantum profile have also been studiedu The constant transverse field H is re- 
placed by H(n) which takes the value +H in the left sector and —H in the right 
sector of the system. Then for S = | and e = 1, the profile e(n) = (0| S x (n) 0) 
can be calculated exactly. In distinction to the monotonous profiles of Figs. ELM 
e(n) is found to show oscillatory behavior in the O phase. For T — + 0, these 
oscillations are suppressed. 



3 Nature of the spectrum 

We have seen that at least qualitatively, the form of the profile m(r) can 
be understood in terms of the subspace M. of the states with a single domain 
wall. While this approximation is a 'priori applicable for h <C 1 ,it is remarkable 
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a) L even b) L odd 

Figure 3: Schematic structure of the low-lying states of "H for S = = 0. 



that the clear separation of the states within M. remains intact even for finite 
values of hu We can thus identify a subspace A4 of low-lying states (with 
M = lim/j^o-M) even though the states in M. are no longer simply those 
states with a single domain wall. Indeed, we find that while the energy gaps 
ofjthe states outside M. are finite within the phases O and F, the states within 
M. show a different behavior. 

We begin with the F phase. Let g(i) = E. L — Eq denote the \ th gap. We 
then findQ the finite-size scaling g{i) ~ L~ e with 9 = 2 for the 2SL + 1 states 
within A4. Furthermore, the scaling amplitudes a(i) = lim£_>oo L 2 g(i) follow 
a simple pattern within M.. Namely {i = 1,2, 3 . . .) 



r(i) = a(»)/a(l) = U{i + 2) = 1, |, 5, 8, y , 16, 



(11) 



We have checked these amplitude ratios within the entire F phase for both 
5=2 an d S = 1. This confirms the earlier observation u that domain walls 
approach their asymptotic behavior algebraically. 

For the O phase, the spectrum is different. We concentratexin S = | and 
h = 0. The basic ^structure of the levels is indicated in Fig. pl3 For L even, 
the states within M. are grouped into doublets with a splitting of order L~ 3 
and the gaps between the doublets are of order L~ 2 . Thus the lowest gap 
scales with an exponent 9 — 3, while 9 = 2 for all other gaps. For L odd, all 
levels are doubly degenerate with gaps between them of order L~ 2 and 9 = 2 
throughout. For the amplitudes, we find independently of the parity of L the 
following pattern, consistent with the structure eq. (flTl) 



p(i)=o(i)/o(2) = 0,l,l,-,-,5,. 



(12) 



c Eqs. dll[|l^ ) can be analytically reproduced by restricting to the space 



M. 
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along the entire line h = 0. Qualitatively, the same structure persists through 
the entire O phase. 
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